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We introduce a new approach to modified gravity which generalizes the recently proposed hybrid 
metric-Palatini gravity. The gravitational action is taken to depend on a general function of both 
the metric and Palatini curvature scalars. The dynamical equivalence with a non-minimally coupled 
bi-scalar field gravitational theory is proved. The evolution of cosmological solutions is studied using 
dynamical systems techniques. 



m 
o 

PL, 



cr 

i ' 



> 
in 
in 
m 
(N 

(N 
O 
m 



x 



I. INTRODUCTION 

In order to explain the present and initial accelerated 
expansions of the Universe a large variety of modified the- 
ories of gravity has been proposed in recent years. Among 
them, one of the most popular is f(R) modified gravity 
where the gravitational action depends on a general func- 
tion of the curvature scalar R, see [l[ for two reviews. To 
derive the gravitational field equations from those mod- 
ified actions, two approaches are extensively used in the 
literature: the metric and the Palatini variational prin- 
ciples (See @ for recent extensions of the Palatini varia- 
tional method in modify gravity). In the so called metric 
approach one takes the metric as the only dynamical 
variable and considers only variations of the action with 
respect to it. The so called Palatini approach is based on 
the idea of considering the connection defining the Rie- 
mann curvature tensor to be a priori independent of the 
metric. As such, one performs variations of the action 
with respect to the metric and the connection indepen- 
dently. Both approaches have been used extensively to 
build cosmological models, many of which contain an era 
of accelerated expansion. 

It is well-known that f(R) theories are dynamically 
equivalent to Brans-Dicke (BD) theories. In fact, metric 
f(R) gravity has been shown to be dynamically equiv- 
alent to a Brans-Dicke theories with vanishing BD pa- 
rameter while Palatini f(R) gravity presents the same 
equivalence if the BD parameter equals —3/2 (see again 
[l|]). The value —3/2 for the BD parameter is a peculiar 
one since it implies no dynamics for the scalar field in 
BD theories. Consequently, Palatini f(R) gravity has the 
same number of dynamical degrees of freedom as general 
relativity, see Q for a very different model which does 
not introduce new dynamical degrees of freedom either. 

More recently, a novel approach to modified gravity 
has been introduced where a Palatini-like f(1Z) term is 
added to the metric Einstein- Hilbert action In this 
context cosmological and astrophysical applications to- 
gether with wormholes geometries have been studied in 
Q, where it has also been shown that viable acceler- 
ating cosmological solutions are allowed by some specific 



models. The theory is dynamically equivalent to a scalar- 
tensor theory with non-minimal coupling to gravity given 
by (1 + 4>)R where R is the metric curvature scalar and 
4> is the scalar field. This is done in strict analogy with 
Palatini f(R) gravity and indeed the BD parameter for 
this theory is still —3/2. However, because of the differ- 
ent coupling to R than in BD theories, in this theory the 
scalar field is dynamical and represents a new dynamical 
degree of freedom. 

In the present paper we analyze a natural extension 
of the theory introduced in We introduce a general 
function which depends on both the metric and Palatini 
curvature scalars. We show that this new generalization 
can be considered as dynamically equivalent to a gravi- 
tational theory with two scalar fields. Only one of these 
scalar fields is non-minimally coupled to R and in gen- 
eral an interaction between the two appears in the action. 
The cosmological features of the theory are then studied 
using dynamical system techniques. 

Let us define k 2 = 87rG/c 4 and start from the action 
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which generalizes the so-called hybrid metric-Palatini ac- 
tion of [4]. In action (flj R is the Ricci curvature scalar 
formed with the Levi-Civita connection, 



(2) 



while 1Z is the curvature scalar of an independent tor- 
sionless connection in analogy with the Palatini ap- 
proach. The variation of the action ([1]) with respect to 
the independent connection leads to 



= o. 
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whose solution is a Levi-Civita connection in terms of the 

1 

2 

The variation with respect to the metric yields 
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where also the matter action has been considered in the 
variation. Because of dU, TZ^ can be related to R^ v 
and terms involving (derivatives of) df /dTZ as in Palatini 
f(R). However the trace of (O now relates TZ and its 
derivatives to T and <? M „, meaning that it is not possible 
in general to solve for TZ. However, this can be avoided 
by requiring d 2 f / dRdTZ = 0. In this particular case the 
trace of ([S]) becomes an algebraic equation in 1Z which 
can be solved for T and (derivatives of) g M „. Note that 
in the following we will keep the function / completely 
arbitrary as it turns out that all results we will obtain 
hold for any (sufficiently smooth) function /. 



II. DYNAMICALLY EQUIVALENT ACTIONS 
AND CONFORMAL TRANSFORMATIONS 



Let us start by considering the action 



S = 



f(<*,P) 
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+ dp [ 
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where a and /3 are two scalar fields. Variation with re- 
spect to a and /3 gives the system 

d 2 f d 2 f 

^(*- a > + <d^-0 =0 ' (7) 



da 
d 2 f 



d 2 f 



d a d ^ R - a) + W {1Z - P)= ^ (8) 

whose only solution is given by a — R and /3 = 1Z, pro- 
vided that 
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This condition follows simply from requiring that this 
matrix type equation is non-degenerate. It is interesting 
to note that the matrix involved is in fact the Hessian 
of /. Since our theory is based on an action principle, 
a non-degenerate Hessian in this context means nothing 
but that solutions of the field equations derived from the 
action are indeed stationary points of the action. 

It is now clear that substituting this solution back into 
action ^ produce immediately action (fT]). The two ac- 
tions are thus dynamically equivalent. Constraint (|9]) 
excludes from our analysis the cases when the function / 
is linear in either a (i.e. R) or f3 (i.e. 1Z). However the 
first case is nothing but the hybrid metric-Palatini theory 
studied in 0, [f| , while the second is equivalent to usual 
metric f(R) theories. Moreover constraint © excludes 
also some particular models such as / = exp(i? + TZ) or 
/ = V 'R7Z. We have thus to reduce the results of this 
section to the models satisfying ([§]). 



Let us define two new scalar fields as 

df(a,f3) , t df(a,(3) 



X 



da 



and £ 



dp 



(10) 



The minus sign in the definition of £ is required in order 
not to allow for a negative kinetic energy of the field. 
Action © can be rewritten as 



(11) 



where the interaction potential is defined as 

V( X ,0 = -/(a(x),0(O) + xa( X ) - £0(0 • (12) 

Due to solution Q (which can also be obtained varying 
action (fTTj) with respect to T^) we can expand TZ and 
find (up to boundary terms) 

1 



S = 



(x-0i?-^(d£) 2 -nx,£) 



(13) 



We can shift x by £ defining a new scalar field as 
X — £• In this way the action becomes 



S = 



(14) 



It is possible to think of action (Tl4"|) as a Brans-Dickc 
theory with vanishing BD parameter and a potential in- 
teracting with another minimally coupled scalar field. 

At this point we can perform a conformal transforma- 
tion in order to switch from the Jordan to the Einstein 
frame. The transformation 



fj,u 



(15) 



allows then, up to surface terms, to rewrite action (TT4"]) 
as 
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Finally, we redefine the two scalar fields as 
13 In (j) , - 2y/2 



2 K 

and the action becomes 



5* = 
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where the new potential is defined as 



HV3 w{e V2^/V3 jK 2^2 /8 y (ig) 



Action (|T5|) is well-known within the context of the 
so-called Brans-Dicke or two-field inflation @43> wnere 
the scalar field <f> represents the Brans-Dicke field while £ 
denotes the inflaton. Usually these studies start from a 
more general action than (|14[) where also a kinetic term 
for <f> is considered, but do not allow for a coupling be- 
tween the two scalar fields in the potential Q. A more 
general coupling between the scalar <fi and R was consid- 
ered in Q. In the Einstein frame this leads to a general 
function of <j> i n the exponential coupling to the kinetic 
term of £. In general we can address action (TT8l) as a spe- 
cific model of Brans-Dicke inflation with vanishing BD 
parameter. This means that from hybrid metric-Palatini 
gravity we have a natural explanation for introducing 
both the inflaton and the Brans-Dicke scalar fields. We 
refer to [8j for recent developments in the context of two- 
field inflation. 

In the next section we will analyze the general cosmo- 
logical dynamics of the theory in the Einstein frame. 



where a(t) is the scale factor. From the gravitational 
field equations ([20]) we obtain the following cosmological 
equations 



3^ + 3H 2 = *Le-^H 2 + ^ 
a 2 - 2 2 



2 +k 2 W . 
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2H + 3H 2 



+ K 2 W. 



(27) 



The scalar fields equations (f2"3"|) and (IM1) give the follow- 
ing evolution equations 

^ + 3^+ * eVW^ + W^O, (28) 
V6 



'2/3k4> 



Wt = 0. 



(29) 



In what follows we will consider only spatially flat 
(k = 0) cosmological models. In order to recast the cos- 
mological equations (f2"6]) -(f29 f into a dynamical system, 
we will make use of the following dimensionless variable 



III. COSMOLOGICAL DYNAMICS IN THE 
EINSTEIN FRAME 

For the sake of simplicity, from now on we will omit 
the tildas in action (|18p. In other words, in what follows 
we will denote with g^, <f>, £ and W the quantities in the 
Einstein frame. The field equations can be obtained by 
varying action (|18p with respect to the dynamical vari- 
ables g^, 4> and £. The gravitational field equations in 
the Einstein frame are thus given by 



n — K 2 ( tW 
where we define 
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(21) 
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For the moment we assume that every other form of mat- 
ter is negligible in comparison to the two scalar fields. 
The equations for the two scalar fields are given by 
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The definitions of the x and y variables have been ex- 
tensively considered to study the cosmological dynamics 
in both uncoupled and coupled dark energy-dark mat- 
ter models Q. With (|30| the Friedmann constraint (|26|) 
reads 



y 2 = i 



impliying that 



< x 2 



< 1 



(31) 



(32) 



since < s 2 < 1. Moreover because of the positiveness of 
the potential energy we must have y > 0, which implies 
that x and y can only take values within half a unit disc. 

In order to complete the autonomous system of equa- 
tions coming from the cosmological equations (|27[) — (f29]) 
we must specify the potential W . In the following we will 
consider three possible cases for the potential energy and 
will analyse the outcoming phase spaces. 
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(23) 
(24) 



where and W$ are the derivatives of the potential 
with respect to (f> and £ respectively. 

We consider a cosmological FLRW metric 



ds 2 = -dt 2 + a(t) 2 



dr 2 



1 — kr 2 



r 2 dn 2 



(25) 



Model 1: W = W e 



First we will consider the usual quintessence exponen- 
tial potential given by 



W = W e 



-\k,<P/ 



(33) 



where Wq and A are both positive constants. In this 
model the only interaction between (f> and £ is given by 
the kinetic coupling being the potential <fi dependent only. 
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TABLE I. Critical points and their properties for model 1. 



Pnint 


X 


y 


Existence 


w eff 


Acceleration 


O IcLLHlliy 


A- 


-1 





V A 


1 


No 


Saddle 


A+ 


1 





V A 


1 


No 


Unstable if A < 6 
Saddle if A > 6 


B 


6 


V2 


A > v/37- 1 


A-2 
A+2 


No 


Saddle if A = V37 - 1 
Stable spiral if A > \/37 — 1 


A+2 


v"A+2 


C 


A 

6 




A < 6 


\ 2 

1 

18 1 


A < 2^3 


Stable if A < %/37 - 1 
Saddle if V37 - 1 < A < 6 


6 



The scalar field <f> plays the role of the usual quintessence 
dark energy field, while we can consider £ as representing 
dark matter. Both dark matter and dark energy have 
thus a geometrical origin in this model and no particles 
have to be introduced opportunely. 

In terms of the variables (f30|) equation (f27|) becomes 

ff = 3(2/ 2 -l), (34) 

which always gives a scaling solution for a(t) in terms of 
the value of y. If y = 1 we have H = and the universe 
undergoes an exponential expansion, while if y > \/2/3 
the universe undergoes a scaling accelerated expansion. 
From ([Ml) we can read off the effective equation of state 
parameter of the total energy content of the universe as 

w eff = l-2y 2 . (35) 

The autonomous system of equation is in this case two 
dimensional and it is given by equations (|28p and (I29|) as 

x = x 2 ~3xy 2 + i(A + 2)y 2 -1, (36) 

y' = -^y(Xx + 6y 2 -6) , (37) 

where a prime denotes differentiation with respect to 
N = In a. There are up to four critical points for this 
system. The points and their properties are showed in 
Table HI There are three possible qualitative behaviours 
of the phase space depending on the following three range 
for A: < A < x/37 - 1, %/37 - 1 < A < 6 and A > 6. 

From Table U we see that in order to have a stable ac- 
celerated attractor we must have A < 2-^/3, so the most 
interesting solutions will belong to the first range whose 
phase space is showed in Fig. [TJ Points A- and A + are 
respectively a saddle and unstable point and represent 
early time solutions with a stiff fluid effective equation of 
state. Every solution evolves eventually reaching point C 
which always lies on the unit circle and acts as a global 
attractor. If A < 2\/3 point C belongs to the region 
above the dashed/red line where the universe is acceler- 
ating. If instead 2-\/3 < A < V37 — 1 point C will be 
below the dashed/red line and the universe will end in 
a decelerating solution. However we can still have a pe- 
riod of accelerated expansion because, for a wide range 



y 




FIG. 1. Phase space for Model 1 with A = 1. The global 
attractor represents an accelerating solution because it lies in 
the region above the dashed/red line. 



y 




FIG. 2. Phase space for Model 1 with A = 4. The global 
attractor does not represent an accelerating solution because 
it lies in the region below the dashed/red line. 



of initial conditions, the evolution still pass through the 
accelerated region for some time as is showed in Fig. [2] 
If A = 3\/2 the global attractor will represent a matter 
dominated universe with vanishing effective equation of 
state parameter, while for A = 2\/6 the final state will be 
a radiation-like dominated universe which suggests that 
this model could be of interest in early time inflationary 
dynamics. 

The phase space for the range \/37 — 1 < A < 6 has 
been drawn in Fig. |3J Point B is now the global attrac- 
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FIG. 3. Phase space for Model 1 with A = 5.9. The global 
attractor is now point B, while C is a saddle point. 




FIG. 4. Phase space for Model 1 with A = 8. The global 
attractor is point B, while now A + is a saddle point. 



tor and always represents a decelerating solution since 
can only appear below the dashed/red line. Point C is 
now a saddle point which attracts all the trajectories be- 
fore they turn to point B. Again, depending on the initial 
conditions, the universe can still undergo a phase of accel- 
erated expansion since several trajectories pass through 
the accelerated region above the dashed/red line. 

The last range for which the dynamics of Model 1 is 
different is given for A > 6. Its phase space is depicted in 
Fig.0] Point A + is now a saddle point and attracts all the 
trajectories along the y direction. The global attractor is 
still point B and the cosmological evolution can experi- 
ence more than one eras of accelerated expansion before 
ending eventually in the final decelerating solution. 

In conclusion we have seen that in Model 1 the universe 
can undergoes phases of accelerated expansion for all the 
possible range of A. If A < 2V3 the cosmic evolution will 
end in an accelerated state, while for all the other values 
of A it eventually reaches a stable decelerating solution. 



TABLE II. Critical points and their properties for model 2. 



Point 



A- 



Aa 



B- 



B 4 



Co 



D 



6 

A+2 



V2 



6 



6 



Existence 



V A 



V A 



V A 



V A 



A > V37- 1 



A < 



A < 



"off 



A-2 
A+2 



9s - 1 



£-1 

18 



Acceleration 



No 



No 



No 



No 



No 



A < 



A < 2\/3 



B. Model 2: W = a «) A e^"^ 
In this section we will consider the potential given by 

W((t>,0 =a( K A e~ AK * /v ^, (38) 

where a and A are two dimensionless positive parameters. 
This potential allows for a direct coupling between the 
two scalar fields <f> and £. Unfortunately, we cannot re- 
cast the cosmological evolution equations (|27|) -([29 |) into 
a two dimensional dynamical system. However defining 
in addition to (I3U1) the new variable 



H 



H + H 



(39) 



we can obtain a three dimensional system. The new vari- 
able z has been choosen in such a way to mantain the 
phase space compact [lfj. It takes values between and 
1, meaning that the phase space is now represented by a 
half cilinder with radius and heigth equal to one. Equa- 
tions (|27p can still be rewritten as (|3~4"|) implying that at 
every point of the phase space the effective equation of 
state is again given by (|35l> . The accelerated region is now 
the part of the half cilinder corresponding to y > y/2/3 
for all the possible values of z. 

The cosmological equations (|27 |) - ([29| give the follow- 
ing three dimensional dynamical system 



x' = y (A + 2 - 6 x) + x 2 - 1 , 



(40) 



y = y 



3 — 3 y x 

y 2 



z' = 3(y 2 -l) (z-l)z, 



2/A 



(41) 
(42) 
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TABLE III. Critical points and their stability properties for model 2. 



Point 


J_J1 tlCll V CLL UCg 


Stahilitv 


A— 


Q O O 1 \ /O 

o, —2, J + Xj I 


Saddle point 


A+ 


3, 2, 3- A/2 


Saddle point 


JD- 


— 3, —2, 6 + A/ 2 


Saddle point 


B+ 


-3, 2, 3- A/2 


Saddle point 


Co 


3A/(2 + A), (-3 ± V81 + 32A - 4A 2 - A ;i )/(2 + A) 


Saddle point 


Do 


A 2 /12, -3 + A 2 /12, -6 + A/3 + A 2 /12 


Saddle point 


Lh. 


-A 2 /12, -9/2 + A/6 + A 2 /8, -oo 


Stable if A < 2(^82 - l)/3 
'Saddle' if (2^ — 2)/3 < A 



where again a prime denotes a derivative with respect to 
N = In a and we have redefine the parameter a as 

P= A 3 ^ a VA ( ^ )2 /A. (43) 

In the above system of equations, the term with /3 in (pLTj) 
becomes singular as y — > or z — > 1 and one must be 
rather careful when investigating the equations for those 
values. The critical points are given in Table [TT] and, 
according to the value of A, there can be up to seven 
critical points. However, when considering the critical 
points with z = 1, we have assumed that the term pro- 
portional to (3 approaches zero when z 1. This issue is 
very difficult to settle analytically, however, the numeri- 
cal solutions and the resulting phase space confirm that 
this assumption is valid. 

Again the three ranges A < V37 — 1 , \/37 — 1 < A < 6 
and A > 6 gives the three qualitatively different be- 
haviours of the phase space. The four points A± and B± 
always represent universes evolving with a stiff matter 
effective equation of state. However they are expected 
to be relevant only at early time and not to be stable 
solutions. 

The phase space for the first range is shown in Fig. [5] 
Points A + and D act as saddle points attracting the 
early time solutions before these turn towards greater 
values of z. The trajectories always evolve towards point 
D\ which represents the global attractor. A few more 
remarks are required about this point. Firstly, one of 
the eigenvalues approaches — oo as z — > 1. The term re- 

/ \2/A 

sponsible for this is /3 y 1 — x 2 — y 2 ( - ^ z _ z ^ J , which 

we assumed to approach zero. The linear stability ma- 
trix (Jacobian) contains the inverse of this term which 
in turn can yield an eigenvalue which formally is — oo. 
This explains why this points act as the global attrac- 
tor to the system. There is a A range where one of the 
eigenvalues of D\ is positive, this point will still attract 
all trajectories. While there is a direction in which the 
point repels trajectories, the non-linearities of the system 
will move any trajectory away from this exact direction 
and the attractive behaviour in the other directions will 
dominate. This behaviour can be seen quite clearly in 
the phase space plots. 




FIG. 5. Phase space for Model 2 with A = 1 and /3 = 0.01. 
The global attractor is point Di representing an accelerating 
solution whenever A < 2\/3 because lies in the region marked 
by the dashed/red line. 



If A < 2\f^ this characterizes an accelerating scaling 
solution, while if A > 2\/Z the universes undergoes dece- 
laration. Again this represent the cosmologically inter- 
esting case where the global attractor of the phase space 
could represent an accelerating universe. 

The phase space for the second qualitative range v37— 
1 < A < 6 is depicted in Fig. [51 Points A + , Co and 
Dq represent saddle points which attract the early time 
solutions. Point Di still represents the global attractor 
but do not characterizes an accelerating solution being 
always outside the accelerated region. 

We have identified another interesting point, C\ with 
coordinates x = y = J^ +2 1 z = 1 which is not a 

critical point to the dynamical system. However, when 
evaluating equations (^H)) and ((32), we note that both 
right-hand sides vanish. The remaining equation (|41[) at 



1.0 



1.0 



1.0 



1.0 



FIG. 6. Phase space for Model 2 with A = 5.9 and /3 = 0.01. 
The global attractor is still point D\ but now the trajectories 
are first attracted by point C\ which is not a critical point. 



this point is 
8 fX 



BC{\) 



l/A-3/2 



\/A 2 + 2A - 36 



1 - z 



1 - z 



2/A 



2/A 



(44) 



where C(A) is a constant depending on A. We can now 
understand the phase space at this point, see Figure [7] 
If 8 is chosen to be small, then the solution behaves as 
if C\ was a critical point. The smaller the value of 8 the 
better C\ acts as an attracting point. However, when 
z gets sufficiently close to 1, the trajectories will get re- 
pelled from this point eventually. If A > 6 none of the 
critical points is stable, and the solution will keep evolv- 
ing without a determined late time behaviour. As can be 
seen in Figure [3 all trajectories reach the z — > 1 surface 
and then stay there. 



C. Model 3: W = W e 



-\k4>/V6 



+ Matter 



In this final section we reconsider Model 1 with the 
potential (|3"3")l and add a standard matter perfect fluid 
energy-momentum tensor to the gravitational field equa- 
tions flu), 



where 



=pg li * + (p + p)U li U v 



(45) 



(46) 



with the comoving four-velocity of the fluid, p its 
energy density and p its pressure. We will consider a 



FIG. 7. Phase space for Model 2 with A = 8 and = 0.01. 



linear equation of state given by 



p = wp, 



(47) 



and assume that standard matter is covariantly conserved 



VT'f) = . 



(48) 



We will also consider only the physical range < z < 1/3, 
meaning that we cannot have cosmic acceleration from 
matter alone. In this model dark matter is included in 
the standard matter sector while both the scalar fields 
4> and £ act as dark energy. Note that we are adding a 
matter action to (I18|) and thus assuming that the matter 
fields couple with the Einstein frame metric g^ v - This 
procedure is different to adding a matter action directly 
to {1} but has largely been considered in literature with- 
out entering in deep phylosophical issues (see [ll] for a 
discussion). 

The new cosmological field equations derived from 
dH]) , with k = 0, read 

3H 2 = n 2 p + if-e-VV^^ 2 + ^ + ^ W , 

(49) 

V, 
(50) 



2H + 3H 2 



- K *p - !£ e -VVZ-- '/I' ... '2l,\- ,.--2 



e-' Y <t> 2 + K 2 w, 



while the two evolution equations for the scalar fields 
still coincide with (f28|) and (l29l) . These equations can 
be recasted in a three dimensional autonomous system 
of equations defining, in addition to pop , the new adi- 
mensional variable 



z 2 = 



2 

up 
3H 2 



The Friedmann constraint (|49p reduces to 



2 i 2 i 2 

x + y + z 



(51) 



(52) 



8 



TABLE IV. Critical points and their properties for model 3. 



Point 


X 


y 


Z 


Existence 


Weff 


Acceleration 


Stability 


A 

J± — 


-l 


U 


U 


V A 


1 
1 


i\o 


C„JJ1 Q 

oaaciie 


A+ 


1 








V A 


1 


No 


Unstahle if A < fi 
Saddle if A > 6 


/ ) 


6 

A+2 


V2 

■yx+2 


n 

u 


\ > ,fxf 1 
a ^_ v o I r 


A-2 
A+2 


Nrv 

LV O 


oaLioic 


C 


A 

6 


V36-A2 
6 





A < 6 


18 ± 


A < 2^ 


Stable if A < 3s/2(w + 1) 
Saddle otherwise 


D 








1 


V A 


W 


No 


Saddle 


E 


3(ra+l) 
A 


A 


x /A 2 -18(tu+l) 
A 


A > 3y/2(w + l) 


W 


No 


Stable if 3^/2{w + 1) < A 
Saddle otherwise 




FIG. 8. Phase space for Model 3 with A = 1 and w = 0. The 
global attractor is point C which represents an accelerating 
solution when A < 2\/3, while point D is a matter dominated 
saddle point. 



implying that the phase space is now the quarter of a 
unit sphere because, thanks to ([521) and the positiveness 
of p and W, we must have y > 0, z > and 

< x 2 + y 2 + z 2 < 1 . (53) 

Equation (|50|) becomes 

§- 2 =\[-2 + 2y 2 + {l- W )z 2 ] , (54) 

from which we can read off the new effective equation of 
state parameter 

w eff = l-2y 2 + (w-l)z 2 . (55) 
The 3D dynamical system is given by equations (|48l) 



Weff 




-1.0 

1 2 3 4 5 6 



FIG. 9. Evolution of the effective equation of state parame- 
ter for the dashed/red trajectories in Fig. [8] After reaching 
a long lasting matter dominated evolution in according to 
observations, the universe eventually ends in an accelerated 
expansion representing the final cosmological stage. 

([S"0]l as 

x' = -^x[2y 2 -(w-l)z 2 ] 

+ ^(X + 2)y 2 + x 2 +z 2 -l, (56) 

y > = -ly[-3( w -l)z 2 + \x + 6y 2 -6\ , (57) 

z' = ^z[(w-l)(z 2 -l)-2y 2 ] , (58) 

and the critical points together with their properties are 
listed in Table IIVI If we compare this with Table U we 
see that we have now two more critical points {D and 
E) corresponding to a universe evolving in accordance 
with the matter equation of state parameter. Point D 
corresponds to a universe completely dominated by the 
matter sector with no dark energy affecting the evolution. 
Point E presents instead both matter and dark energy 
but the total outcome on the universe evolution is still 
completely equivalent to a matter dominated universe. 
The other points, belonging to the z — plane, have the 
same properties of Model 1 with point C being the cosmic 



FIG. 10. Phase space for Model 3 with A = 5 and w = 0. The 
global attractor is now point E where the universe evolves 
according to the matter equation of state. 



FIG. 11. Phase space for Model 3 with A = 5.9 and w = 1/3. 
The global attractor is still point E but now the trajectories 
starting from the z = plane are first attracted by saddle 
point B. 



accelerated stable attractor solution for A < 2\/3. We 
have now four qualitative behaviours for the dynamics of 
the phase space depending again on the possible values 
of A: A < 3y/2(w + 1), 3^2(w + 1) < A < VS7 - 1, 
V37 — 1 < A < 6 and A > 6. Note that because < z < 



1/3 we always have 3y / 2(w + 1) < v37— 1. 

The first range is again the more interesting since we 
can have a late time attractor where the universe is ac- 
celerating its expansion. The dynamics of the phase 
space is depicted in Fig. [8] The late time attractor is 
point C which results in an accelerating solution when- 
ever A < 2\/3. Point D represents a saddle point where 
the universe is completely dominated by the matter sec- 
tor and expands according to radiation/dust solutions. It 
is then clear that every trajectory passing nearby point 
D and eventually ending in point C describes a possible 
physical universe. In fact all these solutions will allow the 
universe to undergo the standard radiation and matter 
eras before the transition to the dark energy accelerating 
solution. 

As an example we can look at the dashed/red solution 
in Fig. [5] and see how the effective equation of state pa- 
rameter evolves. This is plotted in Fig. El We see that the 
universe immediately reaches a matter dominated expan- 
sion and keeps this evolution for some time unperturbed. 
Of course if the matter equation of state parameter w 
changed during this period, from radiation to dust in a 
physical situation, also w e tf would change according to 
w. This means that during this period the universe has 
the time to undergoes the standard cosmological eras in 
agreement with the observations. Depending on the ini- 
tial conditions this stage can last for the time needed to 
produce the nucleosyntesis and create the cosmological 
structures. Eventually we will have the transition to the 
accelerated phase which represents the final phase of the 
universe where the effective equation of state parameter 




-1.0 



1.0 



1.0 



FIG. 12. Phase space for Model 3 with A = 8 and w = 1/3. 
The global attractor is point E but now A+ is a saddle point 
attracting the z — trajectories. 



assumes the value A 2 /18— 1. The situation is completely 
equivalent to Quintessence plus dark matter with (j) play- 
ing the role of the dark energy scalar field. In fact the 
trajectories confined to the border of the sphere, such as 
the dashe/red one in Fig. [U must have s = meaning 
that £ does not influence the dynamics of this solutions. 

We can now look at the second range 3y/2(w + 1) < 
A < v37 — 1 whose phase space is drawn in Fig [10] The 
global attractor is now point E where the universe ex- 
pands according to the matter equation of state. De- 
pending on the initial conditions the trajectories can ei- 
ther pass nearby saddle point D, where we also have 
w e ff = w, or pass through the region surrounding the 
point (0, 1, 0) where the universe accelerates its expansion 
and eventually approaches saddle point C before ending 
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to point E. This shows how this model can be useful in 
early inflationary dynamics since we can have an acceler- 
ated period before the universe start to be radiation/dust 
dominated. 

This feature is presented also in the third possible 
range y/37 — 1 < A < 6 as Fig. Illl shows. For this reason 
and because the phase space properties are more evident, 
we chose to show the dynamics in the w = 1/3 case. We 
still have point E as the global attractor, but now, be- 
sides point C, also point B act as saddle point influencing 
trajectories starting from the z — plane. The dynam- 
ics is similar to the second A-range with several solutions 
experiencing an accelerated phase before ending in the 
matter dominated final evolution. 

Finally the phase space for the last range A > 6 is 
presented in Fig. 1121 Point C is now gone and point 
A + plays its role attracting the z = trajectories. The 
dynamics is again similar to the previous ranges with 
point E being the global attractor and solutions having 
a possible era of cosmic accelerated expansion depending 
on the initial conditions. 



IV. CONCLUSIONS 

In this paper we have studied a natural generaliza- 
tion of the so-called hybrid metric-Palatini gravity intro- 
duced in Q. A completely arbitrary function of both 
the metric and Palatini curvature scalars was considered 
as the Lagrangian density in the action. Using dynam- 
ically equivalent actions and conformal transformation 
techniques, we have shown that this new theory can be 
recast into general relativity plus two scalar fields cou- 



pled with each other. Therefore, using this approach one 
arrives naturally at theories where the different matter 
components couple to each other. It should be empha- 
sized that there are no theoretical restrictions when it 
comes to coupling different matter components, all that 
is required by general relativity and its generalisation is 
that the total energy-momentum tensor is conserved. 

We analyzed the possible applications to cosmology 
considering a FLRW universe and employing dynamical 
system methods. Three specific models specified by their 
potentials were studied in detail and in each case a late 
time cosmological accelerated solution has been found. 
Depending on the model parameters these can represent 
global attractor solutions. We also encountered a rather 
peculiar parameter choice (Fig. [7]) where the dynamical 
system has no global attractor and the cosmological so- 
lution would never stop evolving. 

The first model has a potential without coupling the 
two scalar fields and shows several similarities with usual 
quintessence models. The second model considers a di- 
rect coupling between the two scalars in the potential and 
it is characterized by more mathematical complexity. Its 
evolution is easily understood. Finally, in the third and 
most interesting model we add standard (dark) matter 
to the theory and show that the universe can undergo a 
'extended period' of matter domination followed by an 
accelerating dark energy dominated era. This would in 
principle allow for structure formation in this model. It 
would be interesting to study such models in more details 
studying not only the background evolution but also the 
evolution of perturbation on this background and in par- 
ticular structure formation. This would eventually allow 
us to compare such models with experimental data. 
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